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The fundamental resonances of near-extremal Kerr black holes due to massive scalar perturba- 
tions are derived analytically. We show that there exists a critical mass parameter, /ic, below which 
increasing the mass ^ of the field increases the oscillation frequency 5R(a;) of the resonance. On the 
other hand, above the critical field mass increasing the mass increases the damping rate of 
the mode. We confirm our analytical results by numerical computations. 

O ! I. INTRODUCTION 

(N 

The uniqueness theorems flUsj imply that the metric outside a newly born black hole should relax into a Kerr- 
Newman spacetime, characterized solely by the black-hole mass, charge, and angular momentum. The relaxation 
phase in the dynamics of perturbed black holes is characterized by 'quasinormal ringing', damped oscillations with a 
On ' discrete spectrumfsee e.g. [J, [sj for detailed reviews). These characteristic oscillations are then followed by late-time 
\ decaying tails [1, Q ■ 

■ The black hole quasinormal modes (QNMs) correspond to solutions of the perturbations equations (the Teukolsky 
^ ' master equation [Sj) with the physical boundary conditions of purely outgoing waves at spatial infinity and purely 
ingoing waves crossing the event horizon [9] . Such boundary conditions single out a discrete set of black- hole resonances 
^ ■ {uJn} (assuming a time dependence of the form e~"^*). In analogy with standard scattering theory, the QNMs can be 
f^, regarded as the scattering resonances of the black-hole spacetime. They thus correspond to poles of the transmission 
and reflection amplitudes of a standard scattering problem in a black-hole spacetime. 
I Quasinormal resonances are expected to be excited by a variety of astrophysical processes involving black holes. 

J> ' Being the characteristic sound of the black hole itself, these free oscillations are of great importance from the theoretical 
\ llll and astrophysical point of view [1, i]. They allow a direct way of identifying the spacetime parameters, 

OO . especially the mass and angular momentum of the black hole. This has motivated a flurry of research during the last 
four decades aiming to compute the resonance spectrum of various types of black holes 0, . 
, It is worth nothing that in most cases of physical interest, the black-hole QNMs must be computed numerically 
On ' by solving the perturbations equations supplemented by the appropriate physical boundary conditions. However, it 
, has been shown p^ - [l^ that the spectrum of quasinormal frequencies can be studied analytically in the near-extremal 
] ' limit a — > Af, where M and a are the mass and angular momentum per unit mass of the black hole, respectively. 
. . \ The dynamics of scalar test fields in black-hole spacetimes is primarily of theoretical interest- it usually serves as 
I> ■ a toy model for the analysis of gravitational black- hole perturbations. However, as pointed out in [l^, the possible 
'■ existence of boson stars could make scalar QNMs observationally relevant. Boson stars are assumed to be made up 
''jjj , of self-gravitating massive scalar fields |Ty]. If a boson star becomes unstable and collapses to form a black hole, 
' it is expected to radiate scalar waves (along with gravitational waves) in the appropriate QNMs frequencies. 

Former numerical investigations of massive QNMs (see e.g., [l^lTj]) have found that increasing the mass ^ of the 
field increases the oscillation frequency 5ft(aj) of the mode. Below we shall provide an analytical explanation for this 
phenomena. Furthermore, we shall show that there exists a critical mass parameter, /ic, above which increasing the 
mass jjL of the field actually increases the damping rate 3(w) of the mode. 



II. DESCRIPTION OF THE SYSTEM 



The physical system we consider consists of a massive scalar field coupled to a rotating Kerr black hole. The 
dynamics of a scalar field of mass /i in the Kerr spacetime [3l is governed by the Klein- Gordon equation 

(V'V, - Ai')* = . (1) 

(It is worth emphasizing that ^ stands for A4G/hc, where A4 is the mass of the scalar field. We use units in which 
G = a = h = 1.) One may decompose the field as 

^i,n{t, r, e, cf>) = e™'^5z„(0; a^)i?,™(r; aw)e-*'^* , (2) 
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where {t, r, 9, (p) are the Boyer-Lindquist coordinates flS*], oj is the (conserved) frequency of the mode, I is the spheroidal 
harmonic index, and m is the azimuthal harmonic index with —l<m<l. (We shall henceforth omit the indices I 
and m for brevity.) With the decomposition ( EI, R and S obey radial and angular equations both of confluent Heun 
type coupled by a separation constant K{aLij) [l9|-|23|. The sign of ujj determines whether the solution is decaying 
(w/ < 0) or growing (cj/ > 0) in time. 

The angular functions S{9; auj) are the spheroidal harmonics which are solutions of the angular equation [sl. [20l. [23. 

119 

The angular functions are required to be regular at the poles — Q and 9 ~ n. These boundary conditions pick out 
a discrete set of eigenvalues {Kim} labeled by the integers I and —l<m<l. For a^(w^ — /x^) ^ m? one can treat 
a^(u;^ — /J?) cos^ 9 in Eq. ([3]) as a perturbation term on the generalized Legendre equation and obtain the perturbation 
expansion [2^ 



C50 

~ a^uj^ - /?) + + Cka^\^^ " f^^f (4) 

k=l 



for the separation constants Ki^- The expansion cocfhcicnts {cfc} are given in Ref. [2j 
The radial Teukolsky equation is given by 8, 25, 26] 



a4 f A^^) + + A[2maiJ -K- i?{r^ + a^)]] i? = , (5) 
dr\drJi J 

2Afr + a? and _ff = (r^ + a^)tj — ma. The zeroes of A, r± — M ± {JVP — c?)^!"^ , are the black hole 
(event and inner) horizons. 

We are interested in solutions of the radial equation ([5]) with the physical boundary conditions of purely outgoing 
waves at spatial infinity and purely ingoing waves at the black-hole horizon (as measured by a comoving observer) 



171. That is, 



\er''^'^"^^''^y as r->r+ (y ^ -oo) , 

where the "tortoise" radial coordinate y is defined by dy = [(r^ + a^)/A](ir. These boundary conditions single out 
a discrete set of resonances {w„} which correspond to the quasinormal resonances of the massive field (We 
note that, in addition to the QNMs resonances, the massive field is also characterized by a spectrum of bound states 
[13, [13: HI] which tend to zero at spatial infinity.) 

III. THE QUASINORMAL RESONANCES 

It is convenient to define new dimensionless variables 

_ r - r+ _ r+ - r_ _ uj - 

X = ; T = ; vj = — — ; k = 2ujr+ . (7) 

r+ r+ ^ttTbh 

Here Tbh = 4^(1-'^ 1a^) ^^'^ ^ = ^"^^ ^^"^ temperature and angular velocity of the black hole, respectively. In 

terms of these dimensionless variables the radial equation becomes 

x{x + t)^ + {2x + t)^ + VR = , (8) 
dx'' ax 

where V = H'^/r'^x{x + t) — Kim + 2mauj — p?\r\{x + 1)^ + c?\ and H = r^wx^ + r+fcx + r+wT/2. 

As we shall now show, the spectrum of massive quasinormal resonances can be studied analytically in the double 
limit a — )• M and w — >■ mVL (see [2^ for the massless case). We first consider the radial equation (|5]) in the far region 
X ^ maxjr, M(mr2 — t^^)}. Then Eq. ([5]) is well approximated by 

a;2— + 2x— + yfari? = , (9) 
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where Vfar = (w^ — fj,'^)r^x^ + 2{ujk — fi^r+)r+x + [—Kim + 2mauj + — At^(r^ + a^)]. A solution of Eq. ([9]) that 
satisfies the boundary condition ([5]) can be expressed in terms of the confluent hypergeometric functions M{a, b, z) 

R = Ci(2Vw2 - /^V+)^+^'^a;-5+«'5e-V"'-f''-+^Af (i +i5 + iK,l + 2iS, 2iy/uj^ - y?r+x) + C-iiS -5) , (10) 



where Ci and C2 are constants. Here 



and 



(11) 



k + 2mauj — Ki, 



(12) 



The notation (b — > —S) means "replace b by —b in the preceding term." 

We next consider the near horizon region a; <C 1. The radial equation is given by Eq. (|5]) with V — > Vncar = 
— i^jm + 2TOaaj — /i^(r^ +a^) -V [kx -V vor 12')^ j x[x + t') . The physical solution obeying the ingoing boundary conditions 
at the horizon is given by [23, [11, Isoj 



R = x-^' 



■(7-0 



2F\{— + ib — ik, — ~ ib ~ ik; 1 — iw; —x/t) , 



(13) 



where 2^1(0, b] c; z) is the hypergeometric function. 

The solutions ([TOl) and ([T3|) can be matched in the overlap region maxjr, M(toJ7 — w)} <C a; <C 1. It is worth 
emphasizing that in order to have a non-trivial overlap region we must restrict our analytical solution to the regime 
of rapidly rotating (near- extremal) black holes. In p articular, the condition t <C 1 is satisfied in the near-extremal 
limit. The a; < 1 limit of Eq. ^ yields ^ [H, [sq 



R Ci{2iy/uj^ - fi^r+)^+''^x-^+'^ + C2{b -> -b) . 
The a; > T limit of Eq. dTS]) yields [H, HI 



T{2ib)T{\-ivj) 



x-^2 + ''^ + {b ^ -b) 



r(i +ib - ifc)r(i +ib -ivj + ik) 
By matching the two solutions in the overlap region one finds 



i —iS—i-uj l2 



V{2ib)T{l - iw) 



r(i + ib- ik)T{\ +ib-ivj + ik) 



{2i^i 



up- — /i^r+) 



-i(5 



(14) 



(15) 



(16) 



V{-2i5)Y\\-iw) 



r(i — ib — ik)Y{^ — ib ~ iw -|- ifc) 
Approximating Eq. for x ^ 00 one gets [li, IH, [30l| 



— (2zv/a;2 -^j2^+)-5+'' 



r(i + 2»^) ^_ 

r(i + i(5 - zk) 



C2((5^ -5) 



Ci(2Vw'-m'^+) 



^(-1) 



— i — — i 



+ C2(5 ^ -^) 



(17) 



(18) 



A free oscillations of the field (a quasinormal resonance) is characterized by a purely outgoing wave at spatial infinity. 

Thus, the coefficient of the exponent in Eq. (|18p should vanish, see Eq. ([6|). Taking cognizance of Eqs. 

([TSll - pSp . one finds the resonance condition for the quasinormal modes of the massive field: 



r{2ib)T{l + 2ib){-2iTyUJ^^JPr+) 



-iS 



T{-2ib)T{l - 2ib){~2iT./^J^^^^r+y^ 



r(i -\- ib — in)T{j + ib — ik)T{j + ib — iw + ik) r(i — ib — ZK)r(i — ib ~ ifc)r(i — ib — iw + ik) 



(19) 
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The resonance condition (|T9)) can be solved analytically in the regime r ^ 1 with uj ~ mfJ. We first write it in the 
form 

1 



r(i ^ i5 — ivj + ik) 



V X (^2iTv/a;2_^2^^)-2»5 ^ (20) 



where 2? = [r(2i,5)]2r(i - i(5 - in)T{\ - i6 - ik)/[r{-2iS)]^r{^ +iS~ iK)r{^ + iS - ifc)r(i + i6 -im + ik). We note 
that T> has a well defined limit as a — > M and u — > mil. 

In the limit lu — >■ mfi, where ui is almost purely real, one finds from Eq. ([T^ that is also almost purely real. If 
S is almost purely real and larger than ~ 1, then one has (— i)"^"' = e^~*5 ^'^"^^''^ — e~'^^ ^1. US is almost purely 
imaginary with a positive imaginary part, then one has ^ in the near-extremal t ^ limit. In both cases 

one therefore finds e = (— 2zT-\/a;2 — ^'^r+)~'^^^ ^ 1 on the r.h.s of Eq. (^0)) . 

Thus, a consistent solution of the resonance condition (|20p may be obtained if l/r(i - iS - iw + ik) = 0(e) [3l|. 
Suppose 

i -i(5-in7 + «fc = -n + ?;e + 0(e2) , (21) 

where n > is a non-negative integer and 77 is a constant to be determined below. Then one has 

r(i -iS-iw + ik) ~ T{-n + i]e) ~ {-n)~^T{-n + 1 + rje) ~ • • • ~ [(-l)"n!]"ir(r;e) , (22) 

where we have used the relation T{z -I- 1) = zT{z) [^J]. Next, using the series expansion l/r(z) — X^fc^i '^fe-^*^ with 
ci = 1 [see Eq. (6.1.34) of [H]], one obtains 

l/r(l -i5~iw + ik) = (-l)"n!r/e + ©(e^) . (23) 

Substituting (gg into ^ one finds = ^/[(-l)"?!!]. 

Finally, substituting -nj = (w — mVl) /^ttTbh and A; = 2wr+ ^ m + 0{MTbh) [the last equality holds for uj = 
mn, + 0{Tbh)] into Eq. (PT|) . one obtains a simple formula for the quasinormal resonances of the massive field: 

oj = mil + 27rTBH[m -5~ + i)] + 0{MT^h, ^Tbh) ■ (24) 



IV. NUMERICAL CONFIRMATION 



We shall now verify the validity of the analytically derived formula (IM)) for the massive resonances. The black-hole 
quasinormal frequencies can be computed using standard numerical techniques, see [l3| for details. We present here 
resuhs for the case I ^ m = 2. Substituting in Eq. (g]) a -J> M and Muj = 1-1- 0[MTbh) [see Eq. ([24]) for a -J> M], 
one finds K22 — 6| — jM'^fi^, where we have used the expansion coefficients ci = i and C2 = from [24]]. Next, 

substituting this value of K22 into ((T^ , one obtains 

25 1 

^22 = ^ - lyM'/i' + OiMTBH) . (25) 

One therefore finds that S22 is real for fi < fic — v^25/32Af~^ and imaginary for larger values of the field mass. (It is 
worth emphasizing that in full units M/i stands for the dimensionless ratio GMAi/hc = M Ai / Mp^^^^^.) 

Taking cognizance of Eqs. (f24l) - (|25|) . one finds that for small mass values < /ic, where S22 is real), increasing the 
mass /Lt of the field increases the oscillation frequency ^(ui) of the resonance. On the other hand, for ^ > jic (where 
822 is imaginary) increasing the mass /.t of the field increases the damping rate ^(w) of the mode. 

In Table U we present a c omp arison between the analytically derived massive resonances, Eq. (j24p . and the numeri- 
cally computed frequencies [i7|. We find an almost perfect agreement between the two. Table HIl demonstrates the fact 
that the agreement between the numerical data and the analytical formula (1^^ is quite good already at a/M = 0.9. 
This is quite surprising since the assumption r ^ 1 breaks down for this value of the rotation parameter. 
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5Ra;ana/5Ra;num 




0.0 


1.003 


0.983 


0.1 


1.003 


0.983 


0.2 


1.003 


0.983 


0.3 


1.004 


0.983 



TABLE I: Massive scalar quasinormal resonances of a near-extremal Kerr black hole with a/M = 0.995. The data shown is for 
the mode Z = m = 2, see also [TJ- We display the ratio between the analytically derived frequencies, Wana, and the numerically 
computed values, tJnum. The numerically computed frequencies of the massive field agree with the analytical formula 1)24^ to 
within 2%. 



a/M 


3?tt^ana/3?tt'num 




0.9 


1.007 


1.095 


0.95 


1.010 


1.053 


0.99 


1.005 


0.995 


0.995 


1.003 


0.983 



TABLE IL Massive scalar quasinormal resonances of a near-extremal Kerr black hole. The data shown is for the mode I = m = 2 
with Mfi = 0.1, see also [IJ. We display the ratio between the analytically derived frequencies, ojana, and the numerically 
computed values, Wnum- The agreement between the numerical data and the analytical formula (|24p is quite good already at 
a/M = 0.9. 



V. SUMMARY 



In summary, we have studied analytically the quasinormal mode spectrum of massive fields in the spacetime of 
near-extremal rotating black holes. It was shown that the fundamental resonances can be expressed in terms of the 
black-hole physical parameters: the temperature Tbh and the angular velocity fl. Furthermore, we have shown that 
there exists a critical mass parameter, fJ,c{l, m), below which increasing the mass /i of the field increases the oscillation 
frequency 3?(a;) of the resonance. On the other hand, above the critical field mass increasing the mass of the field 
increases the damping rate of the mode. We confirmed our analytical results by numerical computations. 
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